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Abstract 

We propose a new approach to four-dimensional Planckian-energy scattering in which the 
phase of the 5-matrix is written — to leading order in h and to all orders in R/b = Gs/J — 
in terms of the surface term of the gravitational action and of a boundary term for the 
colliding quanta. The proposal is checked at leading order in R/b and also against some 
already known examples of scattering in strong external gravitational fields. 
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1 Introduction and Outline 



Planckian-energy collisions represent a problem ideally located between the 
interests of particle physics and those of general relativity. At those energies, large 
gravitational fields are both generated and felt by the colliding particles, which can 
then be used as sources and probes for classical and quantum gravity effects. 

By changing the impact parameter b of the collision (or its total angular mo- 
mentum, J = bE) one can explore large, intermediate or even short distances, thus 
progressively increasing the sensitivity to quantum effects. 

This is particularly noticeable when the problem is studied within the context 
of string theory. String theory possesses a fundamental length parameter of its own 

\ s = V^h. (1.1) 

While General Relativity expectations are recovered at large distances 0, substan- 
tial quantum-string modifications occur at b < \ s provided R = GE < \ s . 

On the other hand, for R > X s , string-size effects appear to be negligible flj. In 
this paper we shall assume to be working in the latter regime, so that one is left 
with just two dimensionless parameters characterizing the collision: 

a G = Gs/h and Gs/J = AR/b. (1.2) 

By definition of Planckian-energy scattering, the parameter ac, the gravitational 
equivalent of the fine-structure constant, will always be taken to be large. This 
condition ensures that the scattering amplitude has a large phase, hence that the 
process is semiclassical. We note that this condition also implies that R is much 
larger than both the Compton wavelength of the colliding particles, A c = h/E, and 
the Planck length, lp = \J GTi. 

By contrast, we can still vary at will the second, crucial parameter, R/b. While 
the process is always semiclassical, very different physics is expected to emerge 
depending on the value of R/b. 

The simplest case occurs at R/b <C 1. Here a leading eikonal result has been 
obtained by a variety of methods |], 0, f|, [5], §>). We note in particular, for later 
comparison, the approach of Verlinde and Verlinde [||], where the problem is reduced 
to estimating a "topological" action, that is, a surface term coming from a reduced 
two-dimensional action. 

Unfortunately, the most interesting semiclassical phenomena that should origi- 
nate from the collision — such as gravitational collapse, black-hole formation and 
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Hawking evaporation — are only expected to occur at values of R/b of 0(1) or 
less. This can be seen in various ways, either within an iS-matrix approach (which 
gives [0, H higher-order corrections as a power series in R/b) or from the General 
Relativity point of view (numerical studies of gravitational collapse of rotating sys- 
tems [U , collisions of black holes at high energy and small impact parameters ||. 10 
etc.). 



It has been pointed out [[II] that the approach of ref. || cannot be extended 
in a straightforward way so as to be able to cope with this interesting regime. So 
far, the most promising way to tackle this difficult problem is, in our opinion, the 
one of ref. ]7|, where, following work by Lipatov (T^j, one tries to describe the 
whole series of corrections in R/b in terms of the classical solutions of an effective 
two-dimensional action (in the transverse coordinates). 

The main criticism that one can raise to such an approach is its insistence in 
separating transverse from longitudinal coordinates even in the regime of large angle 
scattering, which necessarily has to precede that of collapse. Besides, the approach 
is technically quite complicated and, so far, only the next-to-leading correction has 
been computed and found to agree with previous direct calculations. 

In this paper we propose a completely new approach, which remains genuinely 
four-dimensional at all stages and yet yields a "topological" result for the 5-matrix. 
This approach is built upon two main ingredients: 

• The correct treatment of the gravitational action through the inclusion of an 
appropriate surface term; 

• A first quantized path integral approach to quantum-field theory developed 
by Fradkin some thirty years ago. 

It makes it possible to express the semiclassical phase, to all orders in R/b, just 
using surface terms. These come partly from the surface term of the gravitational 
action and partly from a boundary term connected to the external particles. 

The outline of the paper is as follows: In section 2 we recall some basic facts 
about the necessity and form of the surface term to be added to the usual Einstein- 
Hilbert action. We also present some general results on how to express such a 
surface integral for asymptotically flat space-times in terms of Bondi masses. De- 
tails of the calculations pertaining to this section can be found in the Appendix. 
In section 3 we use Fradkin's approach to quantum field theory in order to give 



2 



convenient expressions for the full and amputated scalar propagators in an external 
gravitational field. The eikonal approximation is recovered in the relevant limit, 
but we shall not be restricted to it in the following. In section 4 we discuss scatter- 
ing in an external gravitational field, easily recovering known results. We deal, in 
particular, with the case of an external shock-wave and with that of Schwarzschild's 
metric. In section 5 we combine our previous results to obtain a simple expression 
for the iS-matrix of two colliding particles, first in the case of elastic scattering and 
then with the inclusion of gravitational bremsstrahlung. In both cases the process 
is completely determined by surface terms. For the case of elastic scattering we see 
the external metric approximation emerging at leading order (and failing beyond) 
and we recover in an elegant way the leading eikonal result. Section 6 contains some 
remarks concerning the possibility of extending the method to the case in which 
graviton loops are included. 



2 The gravitational action and the <S-matrix of 
massless particles 

Two related facts distinguish the theory of gravitational interaction in an 
asymptotically flat space-time. First, the Einstein-Hilbert Lagrangian is linearly 
homogeneous in the metric, whereas the usual field-theoretic Lagrangians are es- 
sentially quadratic in the field. Second, the energy of gravitating fields is given 
by a surface contribution at the asymptotically flat infinity, whereas, in the non- 
gravitational field theories, it is the integral of a volume density. These two proper- 
ties have not received proper attention in quantum theory. The reason is probably 
that, in perturbation theory, the gravitational interaction is really not too different 
from the others. It seems, however, that any reasonable attempt at non-perturbative 
quantum gravity should rely on the above basic facts. No such attempt will be un- 
dertaken in the present work but we shall argue that, even in the regime where the 
gravitational field is treated classically, the use of the above features of gravity can 
be fruitful. 
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2.1 The gravitational action 

We start by recalling that the Einstein equations with asymptotically flat 
boundary conditions do not follow from the Einstein-Hilbert actionjj] 

& = -i4/ d4 *^ K - <"> 

Indeed, by making a variation of g^ u in a compact domain Q, one finds 
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[ ^xy^gTZ = 
Ju 

jf dV=0 (g^VHg, p - VT^) - (iV" -\gTn) S 9i 



(2.2) 



where V 2 = g af3 V a V/?. The total derivative in ( |2.2| ) can be written in terms of 
contributions from the boundary d O, but when the spatial boundary is pushed to 
infinity with the appropriate asymptotic behaviour of the metric, these contribu- 
tions do not vanish. For the Einstein equations to extremize the action, one should 
add to the Einstein-Hilbert action some functional of the metric to compensate the 
contribution of the total derivative in (512). We denote this functional by Sx[g\- 
The total action is thus of the form 

S = -— ^- / d'x^gn+ S K [g] + Ssource, (2.3) 

where we have included possible matter sources of the gravitational field. The only 
condition defining Sx\g] is that 

SS K [g] = J^x^g-(g^V 2 5g, u - VT^) , (2.4) 

which ensures that the variational equations of the action ( |2.3| ) are of the form 

n^ v --g^TZ = SirGT^, 

q-fiv _ ^ 5 S Source ^ g\ 

Since the action is defined up to an additive constant, we impose also the normal- 
ization condition that Sx[g] vanishes when the metric is flat. 

The derivatives acting on bg^ v in ( [2.4D can be decomposed into normal and 
tangential components with respect to the boundary dQ. Terms with tangential 

1"We use the signature (—,+,+,+), the conventions 1Z utia f) = d v T^ a — 1Z VI1 = Tt^p^ 13 , 
TZ = g^TZ^v , and, in what follows, put c = h = 1 . 
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derivatives can be added and subtracted at will because the action ( |2.3|) is varied 
under the condition 

%,U=0. (2.6) 

For this reason, there are many functionals Sj([g] satisfying the above requirements, 
differing from each other by tangential derivative terms. Addition of such terms 
affects neither the variational principle nor the value of the action fl2.3| ) computed 
with an asymptotically flat metric. The action ( |2.3| ) with Q extended to the whole 
of space-time is thus unique. 

By using the freedom in tangential derivative terms, the functional Sx[g] can 
be put in an elegant form where it is the trace K of the second fundamental 
form Kij on d Q integrated over the boundary. However, we will not be able to use 
this expression, at least directly, because, as discussed below, we have to deal with 
null boundaries. We use instead the simplest expression for Sxlg], the one that can 
be obtained as follows. 

Let the boundary dfl be given by the equation 

dVl : t(x) = 0, (2.7) 

where r(x) is a piecewise smooth function, and let the sign of r(x) be chosen so 
that the direction of growing r leads out of Q. Then the Gauss theorem can be put 
in the form: 

jf d 4 x^ = J d 4 xV=g 5 (r(x)) /"V M r , (2.8) 



which is also valid for null boundaries, (Vr(x)) = 0. Accordingly, eq. (|2.4j) takes 
the form: 

SS K [g] = ^ J d 4 x^ 5(t(x)) V,t (g^WSg aP - g» a V% a p) . (2.9) 

Since the metric enters Sk only near the boundary, and the boundary will be carried 
to the domain where space-time is asymptotically flat, we can introduce an auxiliary 
flat metric g^ v in this domain and expand Sk as 

S K = S K [g + h]= A[g) + B[g, h] + 0[h 2 ] , (2.10) 

where A is the zeroth-order term, B the first-order term, and 0[h 2 ] higher-order 
terms. With the above normalization, 

A[g}=0, (2.11) 
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and the 0[h 2 ] terms fall off at infinity too rapidly to give a finite contribution. The 
result is that, whatever Sk is, it is given by a linear term of its deviation from the 
flat space-time value. But then we already know it, it is: 

Sk = / d A xyf^5(r(x)) V , t r (g^Wh a/3 - g^hj) , (2.12) 

XU/t LJT ''asymptotic 
domain 

where the quantities and operators with tilde refer to the flat metric g^ v , and 

hp, = g^u - 9fiu ■ (2.13) 

The components of g^ u are calculated in a chart covering the asymptotic domain, 
and the components of g^ u are their flat space-time limits. 

The reason for the surface term Q2.12|) to be non vanishing lies in the fact that 



the expansion of the Einstein-Hilbert action in the wave field ( |2.13|) starts with a 



linear term. Hence, the surface term to be subtracted is also linear in the field and, 
finally, the field itself is not decreasing sufficiently fast at large distance for giving 
a vanishing contribution. 

The boundary contribution to the gravitational action has been discussed in the 
literature mainly in the context of Euclidean gravity. We need to calculate it for 
a Lorentzian asymptotically flat space-time. A major requisite of this calculation 
is the determination of the boundary of an asymptotically flat space-time. If one 
thinks of space-time as of a "cylinder" bounded by time-like and space-like hyper- 
surfaces, then, at the space-like portions, the metric will not be asymptotically flat 
at all. The boundary of an asymptotically flat space-time emerges after the Penrose 



conformal transformation, bringing infinity to a finite distance ||14j| . It then consists 
of the spatial infinity (/°), where the space-like geodesies begin and end, the past 
and future time-like infinities (I~ and I + ), where the time-like geodesies begin and 
end, and the past and future null infinities (T~ and X + ), where the null geodesies 
begin and end. In the physical space-time, the 1°, J ± , I ± are understood as infinite 
limits of affine parameters along the respective geodesies. The metric is asymptot- 
ically flat at Z~ , X + and 1°, but not at I~ and I + . An important fact is, however, 
that 1° , I~ and I + are single points, whereas X~ and X + are three-dimensional 
(null) hypersurfaces. Since the surface term in the action is an integral over the 
boundary, we conclude that it is the integral over 

ffl = rul+. (2.14) 
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If space-time has event horizons, the latter also appear as portions of the boundary. 
We shall come back to this point when discussing the iS-matrix. 

More knowledge about the boundary is given by the behavior of the metric in 
its neighborhood. The past and future null infinities are treated similarly. In the 
appendix we shall describe X + following Sachs ||15| , and build the metric near X + by 
considering the congruence u(x) = const, (Vw) 2 = 0, of light rays reaching X + . A 
similar procedure can be followed at X~. In the same appendix, we then calculate 
the boundary term in the action and show that it can be written in the form 



oo roo 



S K = -~( / dnlUu)- I drMJr)) . (2/15) 



oo J — oo 



where M + (u) and M_(t>) are the Bondi masses at X + and X~, the precise defini- 
tion of which is given in the appendix, and the retarded and advanced time are 
normalized by the conditions ( |A.18| ) and ( |A.32j ). 

2.2 The S- matrix of massless fields 

The 5-matrix of the gravitational field coupled to a set of matter fields ip is 
given by the functional integral 

S = J d[ 9tlu ^]e iS , (2.16) 

where S is the full action ( [2. 3D , and the integration is carried out over all fields 
interpolating between the asymptotic fields with fixed operator data. For massless 
conformal invariant fields, these data are at X ± . It is assumed that the measure in 
( |2.16|) includes all gauge-fixing, ghost and 5(0) contributions. 



To leading order in h, the classical path dominates the functional integral, and 
the 5-matrix takes the form 

S tree = exp(iS \ sol , n ) , (2.17) 

where the action is to be calculated on the solution of its variational equations with 
the appropriate boundary conditions. This approximation amounts to a resumma- 
tion of the Feynman diagrams illustrated in fig. 1, as discussed further in section 
5. 

Consider first the case of pure gravity. By virtue of the field equations (|2.5|) , 
the volume density of the Lagrangian in (|2.3|) vanishes, and the action reduces to 
the surface term only: 

S tree = exp(iS K [g] \ soV ^j . (2.18) 
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Then, by using the result Q2.15Q , we obtain: 

(2.19) 



Sfree = exp 



oo 

duM + (u) - I dvM-(v) 

-oo 



which relates the iS-matrix elements directly to the Bondi masses of the classical 
gravitational radiation. 

The vanishing of the Einstein-Hilbert action in ( [2.171 ) is, of course, based on its 
being homogeneous in the metric. However, in general, homogeneity alone is not 
sufficient for the action to vanish. With S((p) homogeneous of degree n in <p, and 
ip infinite-dimensional (e.g. a field), one can only conclude that: 

S ip 

The key point here is that the equations of motion are obtained from variations of 
ip which vanish at the boundary, while the variation of (p needed in the homogeneity 
(Euler) equation for S, being global, does not satisfy this property. With the 
corrected action (|2.3| ) for pure gravity, one arrives at the result (|2.18|) , that is to 



nSiip) = [ dxip(x) — -r-^- + a surface term . (2.20) 

J LP 



the identification of the surface term in ( 2. 20] ) with Sk- We stress once more that 



the reason for the surface term to be non- vanishing in this case lies in the fact that, 
for gravity, the surface term is linear in a field which is not decreasing sufficiently 
fast at large distance. 

If one considers, instead of the iS-matrix, the generating functional for Green's 
functions 

Z[J] = J d[g^, ^]exp (iS + iJ d 4 x^g^J^ , (2.21) 

and chooses the field variables as g^ v or \f—gg llv , the result will again be of the form 
( [2.19|) , with the only difference that the classical equations will be modified by the 



presence of an external source. The total volume density of the Lagrangian in (|2.21 



is again linearly homogeneous in g^, and no additional surface terms appear. With 



other choices of the field variables, the result ( |2.19|) will hold only on shell 



The result ( |2.19|) remains valid if the gravitational field is coupled to non-self- 



interacting conformal invariant matter fields. Indeed, the quadratic actions of these 
fields are reduced to surface terms by their own equations of motion, and, unlike 
the gravitational surface term, they vanish provided one chooses the solutions of 
the source-free equations that vanish at infinity and have finite energies. Since the 
energy-momentum tensor of conformal invariant fields is traceless, the total action 
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is again of the form (|2.18|) although the solution to be used in ( |2.18| ) is, of course, 



different. The absence of self interactions combined with conformal invariance of 
matter fields is a sufficient but not necessary condition for the action to reduce to 
the surface term Sxig]- For example, the action 

- i j d 4 x^g fTV^VW (2.22) 

of a scalar field ip is not conformal invariant but is linearly homogeneous in g^, so 
that (|2.18|) again holds. 



A necessary condition follows from the form of expression Q2.19| ) f\. Note that 
M + (u) has finite limits at both u = — oo and u = +00, and similarly M_(v). The 
limits M + (— 00) and M_(+oo) are both equal to the ADM mass Mq, so that Mq 
cancels in the difference. If matter consists only of massless radiation that comes 
in through X~ and goes out through X + , one has M_(— 00) = and M + (+oo) = 0. 
Generally, for the difference of the integrals in (|2.19|) to be finite, one should at 
least have 

M_(-oo) = M + (+oo) , (2.23) 

which means that the energy carried by time-like sources is conserved separately. 
If this condition does not hold, there should be a non-vanishing volume density of 
the Lagrangian to secure the finiteness of the total action. 

A promising feature of writing the action through the Bondi masses, is that it 
easily allows for the introduction of the machinery and interpretations of classical 



gravity theory, thus making possible the use of solutions || H |10| that go far beyond 
flat-space perturbation theory. The iS-matrix can then be constructed and used, 
hopefully, for studying the semiclassical phenomena discussed in the Introduction. 
In particular, if the relevant classical solution has event horizons, then not all the 
energy entering through Z~ will appear at X + ; a portion of it will fall into the black 
hole. Since the iS-matrix is calculated only between the states defined at Z~ and 
X + , this loss should be felt in (|2.19| ). Alternatively, one may include the horizon as 



a portion of the boundary, which means adding the states defined at the horizon. 



^One of the authors (G.A.V.) is grateful to Don Page for discussing this point. 
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3 Scalar Propagators in Fradkin's Approach 



Fradkin [16] has pioneered an approach to quantum field theory, in which 
integration over a quantum field is avoided by writing the propagators as functional 
integrals of the first quantized theory. Here we use his formulation in the case of 
the propagator for a scalar field in Einstein's theory of gravity and in the presence 
of an electromagnetic field. 



3.1 The Propagator 



The Green function (Feynman propagator) for a scalar field is defined [17 

roc 

G{x,y/g) = (y^H'^x) = i / du(y\ exp [-i(H - ie)u] \x) , 



as 



where, in the case of propagation in an external gravitational field g^ 



2H = -V^V, + m 



-1 



Equation (|3.1|) admits a representation as a path integral as |L8| 



G(x,y/g) = i du 



X{u)=y 



[dX"] [dP^] exp 



dv p ■ x - n 



(3.1) 



(3.2) 



(3.3) 



IX(0)=x 

where the dependence of 7i on P^ is given by replacing the covariant derivative 
V M — > ip^. The correctness of this procedure can also be verified by considering the 
action of a relativistic point in an external metric: 



S = -mJ dXyJ-g^XW. 
This can also be written as 



S 



- I (1A 



-1/2 



g^X") + e l ' 2 m 2 



(3.4) 



(3.5) 



by introducing the auxiliary einbein e. The action (|3.5|) has the advantage of al- 
lowing a straightforward massless limit m — > 0. Eliminating e through its equation 
of motion: 



m 



(3.6) 



one recovers (3.4). 
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According to (3.5), the conjugate momentum is 
so that, upon use of (|3.7|) and (|3.6|) : 



-m 



(3.7) 



(3.f 



The reparametrization invariance of the action (O ) forces the canonical Hamil- 
tonian to vanish so that the full Hamiltonian reduces to the constraint itself, that 
is 



m 



(3.9) 



In the definition ( |3.9| ) we have neglected ordering ambiguities, which give rise 
to a possible additional term proportional to the scalar curvature. The coefficient 
of such a term can be fixed, for example, by imposing conformal invariance for the 
massless, spin-0 wave equation, to obtain 



2H=-V^V u + m 2 --Tl. 

o 



(3.10) 



However, in this paper, the functional integral (|3.3| ) will be used only to the lowest- 
order WKB approximation, where it is unambiguous, and the problems inherent in 
its definition will not be discussed further. 

The propagator is therefore defined by (|3.3|) , which yields, after integrating out 
the momenta, 



G(x,y/g) 
x exp 



du j [ ] /-g(X)dX^ (x - X(0)) *W (y - X{v)) 
/2) jj> (g^X" 



' v m 2 



(3.11) 



The measure in ( |3.11| ) and below must be understood as containing the usual factor 
(27rdz/)~ 2 arising in the integration over [cLP^]. 

The Green function ( |3.11|) can also be generalized in a straightforward fashion 
to include an electromagnetic field by the usual substitution p M — > — eA^ to give 



G(x, y/g, A) — du [y/-g(X)dX*\SV> (x - X(0)) 8^ (y - A») 



x exp 



z/2) / dr (g^X< t X v + 2eA^ 



m 



(3.12) 



The propagators ( |3.11|) and (|3.12|) are connected (in the sense that the vacuum 
diagrams have been already divided out), but not amputated. In order to properly 
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define the iS-matrix a la LSZ, we need the on-shell amputated propagator, in which 
the external legs have been removed: 



G c (p,p'/g,A) = 

lim (p 2 + m 2 )(p' 2 + m 2 ) 



d 4 x d 4 y exp 



ip ■ x — ip ■ y 



(3.13) 

G(x,y/g,A), 



where the Fourier transform is taken with respect to asymptotically Cartesian co- 
ordinates. The prescription we use for obtaining G c directly is the following |. The 
integral over v in ( |3.11|) and (|3.12| ) is replaced by a limit in which the initial and 
final v's go, respectively, to minus and plus infinity. Accordingly, the amputated 
propagator is given by 

1 



G c (p,p'/g,A) 



lim 



lim 

1 2 va — >-oo 



where 



-g(X)dX ft ] exp %A 



A = p- X{yi) - p' ■ X(u f ) + / drL(r 



and 



L(t) 



\ (g^X»X v + 2eA^ 



m 



(3.14) 



(3.15) 



(3.16) 



The prescription ( |3.14| ) gives the same result as the definition ( |3.13| ), as can be 
readily checked in perturbation theory. 

The form ( 3.16Q of the Lagrangian corresponds to the gauge in which e = 1 and 
A = r/m, where r is the particle's proper time. Only affine transformations remain 
as invariances of the action after such a gauge fixing. 

The propagator ( p.!4|) is our starting point. It provides a solution for the quan- 
tum motion of a scalar particle in terms of a functional integral over trajectories. 
Notice that the measure [dX M ] in ( ft. 14 ) includes an integration over the initial 
(X^luj)) and final (X^(uf j) positions. 

At Planckian energies the mass term as well as the electromagnetic interactions 
can be neglected. Furthermore the action A is large and the stationary phase 
approximation should be valid to 0(TT l ). 

t This method, based on a result given in G3], was used within a functional approach in a 
series of papers on scattering amplitudes in QED pQ] . 
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Making A stationary not only with respect to X^{r) inside the integration region 
for t but also, say, with respect to the final position gives: 

X" + T» a (X)X u X a = 0, (3.17) 

X v g lv (vf)=X v {vf)Vr=lfr ( 3 - 18 ) 

At this point the only integration left is an ordinary integration over the initial 
position, everything else being then fixed by the (null) geodesic. As a result the 
bulk term in (|3.15|) vanishes and the phase of the Green function (|3.14j) (phase shift) 
can be computed in terms of the initial and final positions of the particle as: 

G c (p, p'/g) - Lim J d 4 X { exp % [p ■ X { - p' ■ X f (Xi) 

~ Lim J d 4 S exp — i [g • £ + P • A(S) J , (3.19) 

where we have introduced the notations: 

X^Xfa), Xf = X(u f ), 
E = (Xi + X f )/2, A = X f -Xi, 

P=(p + p')/2, q = p'-p, (3.20) 



and we have taken y-g(Xi) = 1 since the asymptotic conditions are defined far 
enough from any gravitation field and we are using Cartesian coordinates there. 
Finally, we have introduced the short-hand notation: 

Lim = lim lim . (3.21) 

„2 p /2_ > _ m 2 u^-oo y _ y 
V j — > CO J 



The final expression for the propagator ( |3.19| ) only depends on Xi and Xf. it 



has thus been reduced to a "boundary" term. The propagator ( |3.19| ) appears to 
take different forms according to which variables are integrated first by means of 
the saddle-point approximation (i.e. by using the geodesic) and which are left to 
the end for a more precise integration method. In (|3.19| ) we have written two 
possible forms for the propagator that we will use in the following. In the second 
one, we have traded the remaining integration over Xi for the one over S, with the 
understanding that A should be expressed in terms of £ and of the external field. 

The final outcome of such a procedure may depend, in general, on the above sep- 
aration of integration variables, the correct choice being dictated by the particular 
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problem at hand. In the following, in order to be on the safe side, we shall always 
imply that our final results, when expressed as integrals, are to be trusted only in 
the saddle-point approximation, where the order of integration becomes immaterial 
and the result is unambiguous. 

3.2 Eikonal Approximation 

As a check of (|3.14|) and of (|3.19| ), we can compute G c in the eikonal limit. This 
corresponds to a saddle-point approximation for the functional integral in which the 
classical trajectory X^{r) is computed to the lowest non-trivial order. In this case, 
the Green function ( [3.14D feels the change in phase along such a trajectory as a 
function of the external field g^i^X) or A^{X). 

Let us consider again the massless case m = in the absence of electromag- 
netism. We write the metric tensor as 

g^v = Vt^u + iKhya, , (3.22) 

where k 2 = 8t[G . 

Solving for the geodesic motion (|3.17f) to lowest order in k yields: 



X»{v) = X»( Vi ) - X u (u i )X ff (u i ) f dr r^(X(r)) (3.23) 

J Vi 



and thus, after use of (|3. 18 



X"{v) = X? + p> - Vi ) - pY T dr [ T dr' T^(X(r')) . (3.24) 

J Vi J Vi 

We can now insert the above perturbative results into ( [3.15| ) to obtain 

A = -q-E + K T V (X(r))p^p v dr + 0(q 2 ) + 0{Kq) . (3.25) 

JVi 

In the argument of we may now use the lowest-order (straight) trajectory: 
X » = ^ + p» [u - (v f + Vi)/2\ , (3.26) 
where 

p» = . (3.27) 

Uf-Ui 
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Inserting these results into ( [3.19| ), the eikonal propagator becomes 

1 



G c (p,p'/h) 



lim 

V f — >oo 



Vf-Vi 



d 4 E 



x exp 



-iq • £ + in J 1 (ll + 



p{r 



Vf + Vi 



(3.28) 



The component of S parallel to p can be easily integrated since the integrand 
depends trivially on it: 



d 4 S exp [—iq ■ E] 



p 



dv exp 



.q-p 
—i v 



d 3 S (. . .) 



(3.29) 



The integration in dv cancels, to leading order in q, the factor l/(vf — fi) in the 
limit defining (|3.14|) and we are left with the integration over the three components 
of £ orthogonal to p. 

After a trivial change of variables we finally obtain: 

" + 0O 



/r+oo 
d 3 b exp -iq ■ b + in i h flv (b + pr) p^p u dr 
J — oo 



where we have introduced 
_ VfXi - UiX } 



(3.30) 



(3.31) 



Vf - Vi 

so that Xq(v) = 6 M + p^v, and the remaining integration is with respect to the 
components of b others than the one in the direction of p. The above result is in 



agreement with the known expression [21 



The electromagnetic case can also be derived along similar lines to give the old 
result 



(3.32) 



G c (p,p/A) ~E d 3 b exp 



iq ■ b + ie A^(b + pr) p^dr 



It is perhaps instructive to give a second derivation of the gravitational eikonal 
starting directly from (|3.14|) and using at first only ( p,17|) at fixed X(vi), X(vf). A 
straightforward calculation then gives: 



A 



-P-A-q-E + 



2(v f - v^ 



(3.33) 



(Vf - Vi) 2 Ju t 



drh 



in' 



£ + A 



T-(Vf + Vi)/2 * 
Vf-Vi 
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The integration over A can be done to lowest order in k by standard perturbative 
techniques. A is replaced in the interaction part of the action by a differential 
operator acting on the remaining integral, which is Gaussian, to obtain: 

Gc(p,p'/g) = Lim J d 4 S exp {-iq ■ S) (3.34) 
—i d d f u f 
X ^ [>/ " ^) 2 dP»dP» 
Terms in ( |3.34J ) obtained by not acting on the exponent can be shown to be sub- 
leading by powers of q/P or of 1/P 2 (uf — Vj) and are thus negligible for small angles. 
One is therefore left with the differential operators acting only on the exponent and 
one remaining ordinary integration variable, which can be chosen to be E. The 
result (|3.30| ) is thus easily recovered. 



-k I drh^E 



exp 



- iP 2( Uf _ v .)/2 



4 Scattering in an External Gravitational Field 

4.1 General Considerations 

Before moving on to the full scattering matrix of the two-body collision, it is 
useful to check our framework in the simpler case of the scattering in an external 
gravitational field. 

We consider two examples for which an exact solution for the geodesies is known: 
the (generalized) Aichelburg-Sexl (AS) shock wave []23| and the Schwarzschild (black 



hole) metric ||24||. Planckian scattering in the Schwarzschild metric has been dis- 



cussed in perturbation theory in [25 



4.2 Scattering by a null shock wave 

The AS metric is an exact solution for the gravitational field produced by a single 
massless particle. Its generalization to an instantaneous light pulse of arbitrary 
energy profile p(y, z) is of the form: 

ds 2 = -dU dV + f(y, z)5(U) dU 2 + dy 2 + dz 2 , (4.1) 

where U = t — x, V = t + x are flat-space null coordinates and the function / of 
the transverse coordinates y, z, the profile function of the shock wave, is related to 
the energy profile p by the Einstein equations: 

Af = -2n 2 p. (4.2) 
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For the special case of a point particle of energy E b (b for beam), this gives the well 
known AS result, for which 



p = E b 5(U) 5 {2 \x T ) 
where xt = (y, z), and 



K 

2^ 



E b \n(y 2 + z 2 



(4.3) 



(4.4) 



The classical trajectory lies in a plane. In a coordinate system in which this is the 
x — y plane it is given, in parametric form, by |26| 



U{y) = 2p u u + U 

V(u) = f' 2 U(u)9(U)/A + f(b)9(U)+V 
x T {v) = {b + f'U(u)B{U)/2,0) , 



(4.5) 



where V(vi) = Vq, U{vi) = Uo and xt{vi) = b are the initial conditions. 

We take z/j = and vj = v to simplify the notation; the propagator can be 
written according to ( |3.19|) as 



G c ( PlP '/AS) 



lim — 

v—>oo v 



d 4 Eexpi {q u E v + q v E u )/2 - q T E T 



+(P U A" + P v A u )/2 - P T A T 

where, from ( |4.5| ), we have: 

1 



(4.6) 



f(b) + -f u (2p u u + U Q ) 



+ V , 



b+-f'[2p u u + Uo 







(4.7) 



and we parametrize the momenta as follows: 

P = + -E-E'costf, -E'smti, 0) 

q = (E'-E, E-E' cos -E' sin 0) . 

Integration over Vq gives a S(q u ) thus enforcing 

E = E'(l + cos$)/2, 



(4.8) 



(4.9) 



17 



so that 

P u = 2E, q v = 2P V = 2(E' - E) and 2P T = q T = —E' sin$ . (4.10) 
Accordingly, the phase becomes 



-q T -b + Ef + 
where U = 2p u v + U . 



{E'-E) + \fE- l - qT f 



(4.11) 



The result (4.11) can also be written as 



- qT .b + Ef + ^[q T -f'E} 2 
by means of the mass-shell condition 
q 2 T = 4E(E' - E) . 



(4.12) 



(4.13) 



The above calculation can be repeated, of course, for any initial transverse vector 
^t(^) — (^1; b 2 ). The Green function is therefore 



G c (p,p'/AS) = Lim 2pU6 ^ f du' I d 2 b exp[-ig T ■ b + iEf(b) 

v Jo J 

x exp iv 1 (qx — f'E)' 
which is dominated by a stationary phase at 
q T = Ef, 



(4.14) 



(4.15) 



corresponding to the known [26] relation between scattering angle and impact pa- 
rameter in the generalized AS metric: 



tantf/2 = -f/2. 
In that saddle approximation we can also use: 

1, 

and therefore obtain for the scattering amplitude 

G c (p,p'/AS) = AE5(q u ) [ d 2 b exp\-iq T ■ b + iEf(b) 



lim — J dz/exp iv' (q^ — f'E) 2 



(4.16) 



(4.17) 



(4.18) 



in agreement with the known result [27 
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4.3 Scattering by a black hole 

The problem of the scattering of a massless scalar particle by a (Schwarzschild) 
black hole is most easily discussed in spherical coordinates. We introduce the Green 
function for the partial waves (E, I, m) and (E', I', m') as 

G(E,l,m,E',l',m') = J d 4 x J d 4 y(l'm'E'\y)(x\lmE)G(x,y) (4.19) 

for which, in standard notation, 

(x\lmE) = Y lm {(f),6)R El {r)exp-iEt 
(l'm'E'\y) = Y i r m ,(<f>,e)R* m ,(r)e W -iE't, (4.20) 

in order to transform to spherical coordinates. 

In this case the propagator for the partial waves can be written as 



G(E,l,m,E',l',m') = lim — - — [ dX* fa) 



v f 

x exp \-iEtfa) + iE'tfa)] 

xY lm (0(^), Y; ml (0(^), efa)) R El (rfa)) R* m , (rfa)) , (4.21) 

where the measure is 

dX^ = dtfa) r 2 fa) drfa) d(f)fa)dcos6fa) (4.22) 

and one integrates over the initial conditions, the final ones being given by the 
dynamics. 



The classical equations of motion are given in parametric form as [ffij 



dr\ 2 i_^_(i_ r Jl 
dv ) r 2 \ r 

d9 b_ 

du r 2 ' 

* 'l-^r\ (4.23) 



dv \ r 

where b = l/E and ro = IMG is the horizon radius. 
We consider a planar motion with m = and 

(j>fa) = 4>fa) , (4.24) 

which implies that the integration over the azimuthal initial angle gives a 5o m >- 
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The integration over the angle 9(vi) can also be performed directly, to yield 
5u> t 21 Pi (cos A6) , (4.25) 
where 



A9 = Ofa) - 0(v f ) . 
Similarly, we can integrate over the initial times to obtain 

J dt{ui) exp[-i(E - E')t(ui) + iE'At] = 5{E - E') exp 



iEAt 



where 

At = t(v f ) - t(vi) . 
In this way we find: 



G(l,E) 



lim 

l/f — * + oo 



1 



r {i>i)c\r{i>i)i exp 



xP,(cosA^) J R ra (r(z/ i )) J R^(r(z/ / )) , 
which, for large I, can be approximated by 



iEAt 



G(l,E)= lim 



i 21 r 

- V; J 



r(ui)dr(ui) 



ia; — * — 00 



x exp 1 

x 



Uf-ViJ r(v f ) 
EAt - p r r(v f ) - p' r r(vi) 



lir 



exp ( UA6 + i— J — exp y—UA9 — i— 



(4.26) 



(4.27) 



(4.28) 



(4.29) 



(4.30) 



where we have taken, respectively, an incoming wave at the initial time and an 
outgoing wave at the final time. In our parametrization, and for Vi and Vf sufficiently 
large, we have r(z/j) ~ r(vf) = L, in such a way that the argument of the exponential 
no longer depends on r(z/j). Therefore, the remaining integration over r{vj) can be 
performed to give 



lr 1 r u J 
/ dr(vi) ~ / dL — > 1 , 

and 2L ~ z// — z/j. The radial momenta are 



(4.31) 



(4.32) 
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The partial wave amplitude (propagator) S(l, E) is thus given by: 
S(l,E) = G(l,E) = exp{2i5 l ) 

with: 

25 t 



E 



— [At - 2L] ± 



7T 



l(A6 + n) + - 



(4.33) 



(4.34) 



Equation ( |4.34| ) is the main result of this section. It expresses the phase shifts 
in terms of At and A8, which are defined in ( |4.28| ) and ( |4.26| ) and can be obtained 
by integrating the classical equations of motion fl4.23|) . 

The free- motion phase shift is just 

25l 0) = ttJ + \ . (4.35) 
Going back to momentum space, the propagator can now be written as 



G(E, ti/Schw) = J2 (2/ + 1) Pi(cos #)G(l, E) 
i 

where d is the scattering angle. 

The sum over I is dominated by a stationary phase at 

±?9= dT 



(4.36) 



(4.37) 



By using the explicit form of 5i, ( |4.34| ), we find the stationary phase condition 



( 4.37 ) in the form [L — > oo): 



±t? = b 



dr 



r 2 y/l — (1 — r$/r) {b 2 /r 2 ) 



— 7T 



(4.38) 



which agrees with the classical equation for the orbit found in p8[ in terms of an 
elliptic integral of the first kind. 

The shift in the time coordinate in (|4.34 ) is particularly interesting. It can be 
written as 

dr 1 



At = Vf — Vi + 2r J 



rmin r-r ^/l - (1 -r / r )(b 2 /r 2 ) ' 



(4.39) 



We see that for r m i n < r$ (corresponding to angular momenta I < Eb, where 
b = 3^/3r /2 is the critical impact parameter below which classical capture takes 
place), At acquires an imaginary part: 



ImAt = 2-KTr 



(4.40) 
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Thus the incoming wave is almost completely absorbed for b < b (recall that, in order 
for our approximations to be valid, we are always at r E » 1), the transmission 
probability being O (e~ 4nr ° E ^J = e~ E l TiI , where 



T, 



H 



1 /47rr 



(4.41) 



is the black hole temperature. 

Finally, the imaginary part (|4.40|) gives an inelastic contribution to partial wave 
cross sections with I < Eb : 



J>+1) 



1 - exp {-E/T H ) 



(4.42) 



The above results are in full agreement with those obtained in the literature 
by a more direct (and involved) computation p9 . In particular, the total inelastic 
cross section can easily be evaluated to give 



7T 
£2 



l=Eb 



£(2* 
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-ixr, 



o • 



(4.43) 



In closing this section we wish to mention that sometime ago Hartle and Hawk- 
ing [[3(J addressed the problem of spontaneous black hole radiation by using an 
approach where complex trajectories dominate the path integral. 



5 Two-Body Collisions 

After having checked our method in the case of scattering in an external field, 
we shall now turn to our real goal, the Planckian-energy collision of two light (mass- 
less) quanta. 

As explained elsewhere [JTL [7|] , the semiclassical approximation should be valid in 
the Planckian energy regime for any finite value of the parameter R/b and provides 
an iS-matrix of the form: 

\nS(E,b) = l —^-F(R/b, logs, log b){l + 0{h)) . (5.1) 

This approximation corresponds to resumming all Feynman graphs which be- 
come (possibly disconnected) tree diagrams once the scalar-particle propagators are 
cut. They are illustrated in figs. 1 and 2 for the case of elastic and inelastic scat- 
tering, respectively. Notice that neither matter loops nor pure graviton loops are 
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allowed in these diagrams. Simple power counting arguments show [I], |7| that these 
loops only contribute to the 0(h) corrections indicated in ( |5.1|) . 

In the functional integral, the semiclassical limit corresponds to a saddle point 
approximation in which the classical solution for and g^ v are inserted in the 
action. The basic idea is thus quite simple: use the homogeneity in g^ v of the bulk 
terms in the gravitational and scalar particle actions to reduce the whole action on 
the classical equations of motion to just surface terms (the gravitational term Sk 
and the term S Pext ). 

We first discuss the case of elastic scattering, reproducing in a simple and elegant 
way the known leading-eikonal result, and then briefly discuss how the method 
retains its validity beyond such approximation. 

5.1 Elastic Scattering 

In the absence of matter loops, the elastic scattering amplitude is given by 

(2vr) V 4 ) ( Pl + p 2 - p[ - p' 2 ) A (p! + p 2 - p[ + p' 2 ) = 

x j[dg flu }exjp(iS g )G 1 c {p 1 ,p[/g)G 2 c {p 2 ,p' 2 /g). (5.2) 

where, for the sake of simplicity, we restrict our attention to the case of two distin- 
guishable scalar particles, thus avoiding complications with crossed or annihilation 
diagrams. 

By means of ( |3.14| ), the right-hand side of eq. ( |5.2| ) can be written as 

hm (^) (Z!) / [J-g( Xl )dXZ] y-g(X a )dX$] [&g»A exp(iS) (5.3) 

for 

S = S g + S 1 + S 2 + S Pext , (5.4) 
where 

S g = -1^2 J V^gK<d 4 x + S K , 

S l = \f Q i &T i g fiV {X i )X!t{T i )X»{T i ), 
S Pext = Pl -X 1 (0)-p' 1 -X 1 (u 1 )+p 2 -X 2 (0)-p' 2 -X 2 (u 2 ), (5.5) 

and Sk is the surface term ( |2.15| ). 
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The classical equations of motion are obtained by a variation with respect to 



K 



x 



K 



£ r dr t tfW (x - ^(r,)) Xf (rOXr(rO (5.6) 

zTo J 



1=1,2 



and with respect to Xf : 

xf + r^(Xi) xfxf = o 



for 



1,2. 



(5.7) 



As we have already discussed in the case of the propagation in an external field, 
the homogeneity of Si and, now, of S g as well, yields an action that reduces, on the 
equations of motion, to a surface term: 



S 



iol'n 



Sk + S, 



Pext ■ 



(5.f 



In the semiclassical limit, we therefore have 

(2tt) 4 <5 (4) {pi + p 2 - v'\ - P 2 ) A (Pi + P2 -> Pi + P 2 ) 

1 /> dX 1 (0)dX 1 (z/ 1 )dX 2 (0)dX 2 (z/ 2 ) exp (iS] 



lim 



pf->0 



Z/lZ/ 2 



(5.9) 



The amplitude Q5.9Q can be rewritten by introducing the variables 

y t = Xi{0)-Xi(vi) 

(X x (0) + X^)) 12 - (X 2 (0) + X 2 (u 2 )) /2 
ip 1 - P ' 1 )/2-(p 2 - P ' 2 )/2 



A 

Q 
Pi 
E 



(5.10) 



after which an integration over £ cancels the overall 5-function on the left-hand 
side of eq. Q5.9Q . We thus obtain 



(5.11) 



A (pi + p 2 -> pi + Pa) 
-1 



lim 

p?^0 



Z/lZ/ 2 



d A exp iq ■ A + iP l ■ y 1 + iP 2 • 2/2 + «<SW 



It is at this stage that (possibly non-perturbative) approximate solutions for the 
classical collision geometries or their generalizations should be used to calculate the 
amplitude ( |5.11 ). In this paper we only carry out a first non-trivial check of the 
method by showing how the known leading eikonal result arises from our surface 
terms. The derivation is then compared with the one using perturbation theory in 
the Newton constant. 
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5.2 Leading Eikonal Approximation 

The leading eikonal approximation corresponds, in our approach, to working 
out the surface terms appearing in ( 5.11] ) to 0(E 2 ). For the external particle 



contribution the calculation is simple since, to this order, one just needs the shifts 
Ui to 0(E). To this purpose we can use the leading order form of the metric which 
consists of the two independent AS shock waves created by the two particles: 

h^ u = + V(2) 



2tt 

-^V,VV v V5(V)ln((y-bf + z 2 ) , (5.12) 

where, in suitable coordinates, the transverse vector b can be identified with the 
transverse components of the four- vector A appearing in ( 5.11|) . 



The calculation of the shifts is now identical to the one used in section 4.2, in the 
case of scattering from an external field. In a generic boosted frame, in which the 
two particles are coming against each other along the x axis, the initial momenta 
are 

pi ~pi = (tfi,^!, 0,0) (5.13) 

and 

p' 2 ~p 2 = (E 2 ,-E 2 ,0,0) (5.14) 
and we readily find (see eq.(fOj): 



P 1 -y 1 =P 2 -y 2 = - f ^^\nb 2 . (5.15) 

We see that the contribution of the external momenta is twice the desired value 
because, unlike in the external metric problem, there are now two identical contribu- 
tions from each colliding particle. If the known result is to emerge, the gravitational 
surface term Sk should contribute as: 

S K = ^^\nb\ (5.16) 

Z7T 

We shall first show how this result arises very simply from the Bondi masses 
and next compare it with the perturbative calculation. An apparent difficulty in 
both approaches is that we need the action to 0(E 2 ) whereas the metric ( |5.12j ) is 
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valid only to to 0{E). Nevertheless, as we shall see, having the metric to 0(E) is 
actually sufficient. 

We recall, from the appendix, the meaning of the Bondi masses. M_(t>), at X~, 
is the energy brought into the system by the (advanced) time v, while M + (u), at 
X + , is the energy present in the system at (retarded) time u. Note the difference in 
notation and meaning: the arguments v and u of the Bondi masses are exactly null 
variables and label, respectively, converging and diverging spherical congruences, 
whereas the V and U coordinates appearing in ( |5.12| ) are only approximately null 
(they are such only in the flat metric) and label plane congruences. 

To lowest order in Ei, E2, the only relevant effect is again the above mentioned 
shift of the geodesies computed this time in terms of the variables u, v. Consider 
first the case of a single non-interacting massless particle. In this case, U remains 
exactly null and the Bondi masses are of the form: 



expressing the fact that a particle of energy E\ comes in through X~ at the instant 
v = and goes out through X + at the instant u = 0. There are no further energy 
fluxes either through X~ or through I + . 

In the case of two particles, there are still no interactions at past null infinity. 
Therefore the Bondi mass M_ is simply 



which means that both particles show up at Z~ at the instant v — 0. 

The time instants of their arrival at X + will now be shifted in comparison with 
(|5.171) . By using the geodesic shifts in U and V and the asymptotic relation between 
the two sets of coordinates, it is easy to show that the first particle will arrive at 
J + at u = Ui and the second at u = u 2 where 



M_(«) = E x 9(v) , M+(u) = E l 6{-u) , (E 2 = 0) 



(5.17) 



M_{v) = E x 6{v) + E 2 6(v) 



(5.18) 



K 



2 E 2 



ln& 2 , 



2tt 
■ 2 E 



Info 2 . 



(5.19) 



The Bondi mass at 1 + is, therefore, 



M. 



(u) = E x 9( Ul -u) + E 2 Q{u 2 - u) . 



(5.20) 
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It should be stressed here that, while ( |5.18|) is an exact expression, ( |5.20|) ac- 
counts only for the lowest-order effect of the time shift and neglects higher-order 
contributions from gravitational radiation and mass renormalization. 

Inserting eqs.( |5.18"|) and (|5.20| ) into the expression Q2.15| ) for Sk' 

1 f+oo r , 

S K = -- dx M+(-x) - M_(x) , (5.21) 

2 J—oo ^ 

we obtain precisely the desired result (|5.16| ). Note that the time shifts (|5.19|) are 
to be computed from the metric only to O(E) since there is already a factor E in 
the Bondi masses. This is the way a metric of O(E) leads to an action of 0(E 2 ). 
With Bondi masses limited to first order in E, one would have for M + an expression 
similar to ( |5. 18| ) and a vanishing result for Sk, in complete agreement with the fact 
that there is no contribution linear in E in the action. 

The perturbative calculation proceeds instead as follows. The surface term Sk 
is transformed back to the volume integral of a total derivative (see ( [2.12| ) after 
using Gauss' theorem) yielding 

Sk = J^L d ' x ^ (^W " V"V"V) > (5-22) 

where the tilde stands for flat space-time. 

For a direct calculation of ( |5.22j ) we would need a solution to second order in E 
since, again, the insertion of the lowest-order solution ( |5.12| ) in ( |5.22j ) yields zero. 
This problem can be avoided, however, by noticing that for a traceless source — such 
as our two massless particles- 



TZ=(V z ~g 



[IV I 







(5.23) 



and therefore Sk can be written as the quadratic part of the Hilbert action plus 
higher-order terms, that is 



Ik 2 Ja 



s K = — & a xJ^ 9 k v \r a g u ^ 2 + r^v 



(5.24) 



The metric ( |5.12|) is now sufficient because ( |5.24|) is quadratic in the solution. 
We thus obtain that 

k 2 Ei E 2 



Sk= 8^I ^'-"f-S V (s m 9 M1 ^) hat 



— " l62 ' 



(5.25) 
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since in this case 



aj3 



0. 



(5.26) 



We have thus obtained, by two conceptually different calculations, the same 
final result for the surface action Sk- Comparison of the two methods confirms the 
validity of our approach. It also shows, in our opinion, that working in terms of 
the Bondi masses is simpler and physically more meaningful. The crucial test will 
come, of course, with non-perturbative computations, something we do not want to 
attempt here. 

The scattering amplitude (|5.11| ) can now be computed. In order for the two 
particles to scatter at all, we have to impose the two conditions 

X-{u 2 ) > X-(0) , X+( Vl ) > X+(0) . (5.27) 

This gives two ^-functions that can be written in terms of the variable A as 



9 (A + - E t v-^ 9 (A" - E 2 v 2 ) . 



(5.28) 



The condition ( 5.28|) enforces the requirement that between v = and v = V\, as 
well as between and u 2 , the two trajectories cross each other. The integration 
over the light-like components of A thus gives 



dA + dA" exp 

exp (±q + EiVi 



l -q + A + + ~g_A" 



exp 



9 (A+ - Eiui) 9 (A - - E 2 is 2 
{{q-E 2 i 



iq+ 



iq. 



E 1 E 2 v x v 2 (5.29) 



canceling the factor \jv\V 2 in ( |5.9|) . We therefore obtain: 

r r ,.2 
A (pi + p 2 -> v'x + P' 2 ) = EiE 2 J d 2 6exp 

in agreement with previous results |], ^|, ||, || 



.K 2 EiE 2 2 

iqT ■ o — % m b 

2tt 



(5.30) 



5.3 Inclusion of Gravitational Bremsstrahlung 

In order to extend the calculations of the previous section to higher orders 
in R/b, one has to deal with the problem of gravitational bremsstrahlung (see 
fig. 2). Equivalently, in the language of particle physics, one encounters infrared 
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divergences in amplitudes with a given number of final gravitons and some Block- 
Nordsiek-type treatment becomes necessary in order to be able to extract finite, 
physically meaningful results 0. 

The general reasoning given in section 2.2 insures that the surface nature of the 
semiclassical 5-matrix remains valid even after allowing for the effects of gravita- 
tional radiation. It is however instructive to see how this works in detail in the 
case of two-body inelastic collisions within a conventional LSZ treatment of the 
iS-matrix. 

According to LSZ the amplitude for two scalar particles (if)) to scatter with 
emission of n gravitons (h) is given (symbolically) by 

A(2 ^2 + n) = (if) if) if) if) hh ■■■ &)| LSZ > (5.31) 

that is by applying the standard LSZ steps (Fourier transform, truncation of prop- 
agators and going on shell) to the 4 + n-point function shown in (|5.31 ). 
These amplitudes are generated by the functional 

A(2 -> 2 , J) = (if) if) if) if) exp(/i J)) | Lgz , (5.32) 

through standard functional differentiation with respect to the external source J. 
In the absence of scalar-field loops, the LSZ operation amounts to replacing pairs 
of scalar fields by the corresponding amputated propagators G c and the source J 
by A^h as (x), where A^ 2 ^ is the quadratic operator appearing in eq. ( |5.25|) . h as (x) 
then becomes the new external source. 

Finally, in order for the produced gravitons to be on-shell, the condition 

A®h as (x) = (5.33) 

has to be imposed at the end of the computation. This allows us to replace h(x) in 
eq. ( |5.32j ) with the full g^x). 

We conclude that a power expansion in h as of the iS-matrix functional 

W(h as ) = J [dflv] exp [iS g + iSbr] G^pup'J g)G 2 c (p 2 ,p' 2 / g) , (5.34) 

where 

S br = J d A xg(x)A^h as (x), (5.35) 

gives the amplitudes for emission of any number of gravitons. Following the proce- 
dure of Section 5.1, we can rewrite ( |5.34j ) as a functional integral over g^, Xf and 
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X% with a modified action S containing, beside the terms in ( |5.4| ), the additional 
bremsstrahlung term Q5.35Q . 

The crucial observation at this point is that the additional term ( |5.35| ) is still 
effectively homogeneous of degree one in g^ v . Therefore, within the stated approx- 
imations, we are still left with the result (15.81) and 



W(h as ) = exp [i (S K + S Pext )\ sol , n \ , (5.36) 

as expected from the general arguments of section (2.2) 

Note that, while the two surface terms occurring in ( |5.36| ) are still the same ones 
as in the elastic case, their actual values now include, implicitly, a dependence on 
gravitational radiation through h as (x). Giving h as (x) amounts to specifying a defi- 
nite final coherent state of gravitons, the closest one can get quantum-mechanically 
to a classical radiation field. For a generic h as (x), one will still be plagued by in- 
frared problems. However, as in the analogue QED problem, a well-defined coherent 
state for which infrared divergences cancel will be chosen by the dynamics and will 
represent quantum mechanically the classical asymptotic gravitational radiation 
field generated in the collision. 

In concluding this section, we want again to remark that the validity of ( |5.8| ) 
does not depend on any approximation except the tree approximation, for both the 
gravitational field and the first-quantized particles. The surface term has one and 
the same form for any interactions of the gravitating fields but its value is different 
for different interactions. 

To what extent the surface nature of gravity persist when graviton loops are 
taken into account is discussed in the next section. 



6 Inclusion of graviton loops 

The extension of the above ideas beyond the classical approximation for the 



metric is based on |32|, |33], [34] which dealt with the quantum fluctuations of the 
conformal mode of the gravitational field. In this section we first give a new, more 
direct, derivation of the result of these works that also shows at which points care 
should be taken when using these results; next, we discuss its relevance for the 
present work. 
Let us write 

g^v(x) = g^(x)(j) 2 (x) , (6.1) 
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singling out the conformal mode of the gravitational field, and impose on g lll ,{x) a 
constraint 



X(9(x))=0, (6.2) 
such that the equation 

X (tf(x)g(x)) = 0, (6.3) 
with g satisfying ( |6.2|) , has only one solution 

n 2 (x) = 1. (6.4) 

We may then reparametrize the ten components of g^ v into nine independent com- 
ponents of g^ V) constrained by condition (|6.2|) , and the conformal mode 4>(x). 

With respect to g^ u , eq. ( |6.2| ) plays the role of a gauge-fixing condition, which 
removes the arbitrariness of local conformal transformations g^ u — > Q 2 g^ u . We shall 
consider only the transformations that become the identity at the asymptotically 
flat infinity; the boundary conditions for the conformal mode are 

0(*)U I+JO = 1 • (6-5) 

The choice of the gauge- fixing function x{g) is limited only by the requirement that 
the continuous matrix 

Q(x,y) = M*)^7T ( 6 - 6 ) 
be invertible. 

Our goal is to rewrite the functional integral ( |2.16| ) in terms of a set of variables 
containing the conformal mode <fr(x) explicitly. However, instead of introducing </> 
and the nine independent components of g^ u , we introduce eleven variables — <fi and 
the ten components of g^ v — and insert the delta-function <5(x(<?)) i n the measure. 
By calculating the relevant Jacobians, we find 

1] dg, u = (det Q) J[ 19 5(x(s))d0 J[ dg, u , (6.7) 



with Q given by fl6.6| ). The measure in ( |2.16| ) contains local factors in g^ u [^TJ which 



after the substitution of ( |6.1| ), will change the total power of <fi in the transformed 
measure. However, in what follows, we shall ignore the local factors in the measure 
altogether by assuming, for instance, dimensional regularization. 
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Next, let us find the form of the action 



in the new variables. We have 



(6.1 



where the quantities and operators with a bar refer to g^ v . The transformation of 
the boundary term Sk\q\ is also non-trivial, and its effect is that the total action 
contains no second derivatives of 6: 



^ J d 4 x^n + S K [g] 
-|/dV=I 



1 



-(V0) 2 + - n(g)<p 2 



so that its variational equation in is of the form 



S 



V 2 



0. 



(6.9) 



(6.10) 



As a result, the functional integral ( |2.16| ) for the case of pure gravity takes the 
form 



S= /[d^][d0](detg)<5(x(^))exp -- d 4 x^g -(V0) 



3? 



6 



+ iS K \g]} ■ 



(6.11) 



The conformal mode can now be integrated out. Because of the boundary condition 
rBT), one must first make a shift 



(6.12) 



i- ,i+ ,i° 



<j>(x) = M x \g) + v( x ) > 

where <j)o(x\g) is a solution of eqs. fl6.10|) , ( |6.5|) , and 

(p{x) _ _, = 0. 
We have 

M x \g) 

where G is the Green function 
G 



1 + -GTZ 

6 



ys+gTl 



(6.13) 



(6.14) 



(6.15) 



whose boundary conditions, apart from the fact that the contraction G TZ vanishes 
at infinity, are to be discussed below. 
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Equation ( 6.11 ) takes the form 



S 



[dg nu ] [dip] (det Q)5(x(g)) exp\ iS conf [g] 



K z 



with 



S C onf[g] 



1 



2k 1 



d 4 aV^ K(g) + 7 H(g)GH(g) 



6 



+ S K [g\ 



(6.16) 



(6.17) 



The crucial point is the integration over the field ip. Its action has a wrong sign, 
which is a consequence of the fact that, with Euclidean signature of the metric, the 
Einstein-Hilbert action, even with the term Sj([g] added, is unbounded from below 
|34]j . It is the conformal mode that makes it unbounded. The Euclidean functional 



integral over cp diverges exponentially and we must rotate the integration contour to 
imaginary ip [34|. In the Lorentzian context, the prescription to be used is different: 
one must take for the propagator of the field ip, instead of the Feynman's Green 
function its complex conjugate: 



V 2 + 

6 



l( 



(6.18) 



One thus obtains 

S = J [dfa] S( X {g)) (detQ) (det( V 2 + \^~ «0) f exp(iS conf [g]) . (6.19) 
The final step is the choice of the gauge-fixing function x{g) satisfying the con- 



dition detQ 7^ 0. Equation ( |6.19|) is valid with any such x(<?)> but let us choose 
X(g) = K{g) . (6.20) 



Then, for Q in eq. ( |6.6| ), one obtains 
Q(x,y) = 3(^-ln) 5(x,y). 



(6.21) 



Under the condition 1Z = implied by the delta-function in (|6.19|) , this guarantees 
that detQ 7^ 0. Furthermore, under the condition 71 = 0, one has from ( |6.17f ): 



Sconf [g] 



S K [g] 



and the functional integral (|6.19|) takes the form 
S = [ [dg uu ] 8{K{g)) (det V 2 ) 1/2 exp(iS K \g] 



(6.22) 



(6.23) 
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with only the surface term in the action. 

The action S con f[g] in (|6.17|) was first obtained in [32] and next discussed in 
[ |33| . 34 1 . It is invariant under local conformal transformations that become the 
identity at the asymptotically flat infinity: 



S C onf[^ 2 g] 



Sconf [g] 



(6.24) 



Indeed, the action S con f[g] is obtained from the action of Einstein's theory by making 
the substitution 



9nv{x) = g^ v {x)4>l(x\g) 



9 = 90„(9) 



Sconf [g] = S[g] 
S[g] = ~ J d'x^n(g) + S K [g\. 



Since the function (j) (x\g) in ( |6.14|) transforms as |32j, ^3[ 
(f) (x\n 2 g) = n^ 1 (x)(j)o(x\g) , 



(6.25) 
(6.26) 

(6.27) 
(6.28) 



the combination on the right-hand side of ( |6.25|) is conformal-invariant, and hence 
any functional of this combination is conformal invariant. The action S con f[g] 
may be regarded as a conformal-invariant part of Einstein's action. By rewriting 



eq. ( 6.17 ) in the form 

%] = S conf [g] + J-j | ^x^KGK, 



(6.29) 



one makes it obvious that Einstein's action is the broken action S con f with the 



gauge-fixing function ( |6.20|) introduced in the action quadratically. Since the action 
S con f[g] does not contain the harmful conformal mode, one may expect that, with 
Euclidean signature of the metric, this action is already positive-definite. This is 
indeed the case. In fact, by using the conformal invariance, the action S con f[g] can 
be brought to the gauge 71 = where 



Sconf [g] = S K [g] 
and it has been proved [^] that 



(6.30) 



Sk [g] 



> 0, sign g = (+ + ++) . 



(6.31) 
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In the derivation above, the final result in eq.( |6.23|) appears to be obtained by 
identically transforming the original functional integral. However, even apart from 
the problem with the propagation of the conformal mode, one reserve remains: the 
bar on the integration variables g^ v in eqs. ( |6.19|) and (|6.23| ) is kept because the 
asymptotic fields for and g^ v need not be one and the same. Equation (|6.19|) 
looks like the result of the quantization of the invariant action S con f obtained by 
the usual rules of gauge theory with S(x) the gauge-breaking term and detQ the 
ghost term. The only apparent difference from the usual rules is that, because the 
action S con f is non-local, there appears an additional invariant measure 

/ 1 \\~ 1/2 

det V 2 + - K - ie)) . (6.32) 



6 

However, gravity theory is not conformal- invariant, and the choice of the conformal 
gauge is relevant for the physical metric. There is no other place in the integral 
( |6.19D except the asymptotic condition where the conformal invariance can be bro- 
ken. On the other hand, there is one and only one choice of the conformal gauge 
with which the metric measures distance. For pure gravity or gravity coupled to 
conformal-invariant matter fields, this choice is 1Z = 0. At the classical level, the 
variational equations of the action S con f[g], supplemented with the equation 1Z = 0, 
are equivalent to the Einstein equation. It is, therefore, plausible that, for g^ v in 
( |6.23| ), the original asymptotic conditions apply and, generally, g^ u in ( |6.23| ) can be 
identified with the operator of the physical metric. The functional integral ( |6.23|) , 
or its generalization to the presence of massless matter, could then be used directly 
for calculating the expectation values of the gravitational-field observables and the 
scattering matrix. 

The relevance of the result fl6.23|) is that, after integrating out the conformal 
mode, the action in the functional integral becomes a pure surface term. Therefore, 
this action can be written again in terms of Bondi masses, those pertaining, this 
time, to the virtual gravitational fields that are summed over in the functional 
integral. By extending the approach advocated in this paper one can try to use 
the above result to account for the contributions to the scattering amplitude from 
non-perturbative virtual field configurations. 
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Appendix 



A Calculation of the boundary term 

Let us introduce on each null surface, where u is constant, two coordinates 
0°, a = 1,2, which label the light rays, and a third coordinate r, which is a pa- 
rameter along the rays. The coordinates <p a take values on a 2-sphere and can be 
chosen so as to ensure orthogonality, that is (Vm, V0 a ) = 0. The parameter r can 
be chosen to be the luminosity distance along the rays by considering the induced 
metric f\ on u = const : 

ds 2 \ u=const = g ab d<j) a d<j) b , (a, b = 1,2), (A.l) 

where 

9ab = 9ab(4>,r,u)\ u=const , (A. 2) 

and requiring that the area of the two-dimensional section u = const, r = const be 



/ 



d 2 <f)^detg ab = Anr 2 . (A.3) 

2 — sphere 

The functions u, r, <p a can serve as local coordinates. Four components of the 

metric in this coordinate frame are already known. Three of them are 

g uu = (Vu) 2 EE (A.4) 

g U(pl = (V«,V^) = 0, (A.5) 

g ue = (Vu, V0 2 ) = O, (A.6) 



and the fourth is fixed by the constraint ( A.3 ). The remaining six components of 
the metric characterize the gravitational field. We introduce six unknown functions 
W, 7, 5 and U a (a = 1, 2), so that 

g ur = (V«, Vr) = i (A.7) 
ST = (Vr) 2 = ^W, (A.8) 
g rr = (Vr, V0 a ) = i C/ a (a = 1,2), (A.9) 



^Thc induced metric on u = const is two-dimensional because, along the rays, the interval is 
zero. 
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where 



? = 6. 



ip < < n, <cp < 2ir 



(A.10) 



and 



(Vr, V0) 



U\ (Vr, Vip) = . 



Next 



(V0 a , V0 b ) = g ab = (g 



'ab ) 



where 



9ab> 



- e 
2 



2T 



+ e 2<5 )d# 2 + (e 



7—5 



sin d 9 d 



sin 2 



and 



det = r 4 sin 2 6 1 , 
in agreement with (|A.3| ). Then the metric takes the form 



Wdu 2 + 2^dudr + g ab {d(t) a - U a du) (d<p b - U b du) 



with g a b given by ( |A.13| ), and its determinant is 
g = det = -* 2 r 4 sin 2 6>. 



(A.H) 



(A.12) 



(A.13) 



(A.14) 



(A.15) 



(A.16) 



As seen from ( |A.16| ), the chart u, r, <p a breaks down when \l/ becomes zero or 
infinity. The reason for that can be understood by considering the null geodesic 
u = const, <p a = const. Let A be an affine parameter along this geodesic, growing 
towards the future. From the geodesic equation one finds 



dr 
dA 



u— const 
(j>=const 



oc (Vu, Vr) 



1 



(A.17) 



from which it follows that, at the point where \l/ becomes zero or infinity, r stops 
being a monotonic parameter along the light ray passing through this point. At this 
point the light ray hits the apparent horizon. If, during the history of an outgoing 
light ray, dr/dX changes its sign an odd number of times or if dr/dX — ► as A — > oo, 
this light ray will never reach the asymptotic domain of infinite luminosity distance, 
where the metric becomes flat. By definition, this light ray will not come to X + . 
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By considering only a portion of the full congruence of null geodesies that come 
to X + , we automatically guarantee that the chart u, r, <j) a covers the asymptotic 
domain near X + . Equation ([A .15 ) gives the general form of an asymptotically flat 
metric near X + , and X + is reached at the limit r — > oo with fixed u and <p a . The 
remaining arbitrariness in the choice of the variable u is fixed by the normalization 
condition 



(Vu, Vr) \ x+ = -1 



Then 



(Vu, Vr) = ^ < 



(A.18) 



(A.19) 



everywhere in the domain covered by the chart u, r, <ft a . The u and <p a act as 
coordinates on X + , which is topologically a product of the time axis by the 2- 
sphere. Because of ( |A.18| ), the retarded time u grows along X + towards the future 
and coincides with the proper time of an observer at rest at large and constant r. 
The metric at X + behaves as follows (r — > oo) 



W 



1 + 5 



7 — 5 



+ o(- 



X+ 



efface) 



i+ 



+ o{- 2l . 



x+ 



2A^K y, 9) (I 



+ ' 



(A.20) 
(A.21) 
(A.22) 
(A.23) 
(A.24) 
(A.25) 



where + , C+, C%, J\ff,J\f+ are finite functions of u, <p, 6. 
According to ( |A.8| ), 



(Vr) 



2GM + (u,ip, 



i+ 



+ O 



(A.26) 
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and A4+ averaged over the two-sphere gives 
1 



MAu) 



r 2n 



d#sin# M+(u,<p,9) 



(A.27) 



which is the Bondi mass at X + , that is, the energy that remains in the system by 
the instant u of retarded time. Its limit at u —>■ — oo, 



oo 



(A.28) 



is the ADM mass which is the full conserved energy stored in space-time and 
measured at spatial infinity 1° . The difference 

dM, \ 



MAu) 



du 



du 



(A.29) 



is the energy radiated away through X + by the instant u of retarded time. This is 
a sum of the energy carried away by the gravitational waves and the energy radi- 
ated by massless sources. The energy conservation law following from the Einstein 
equations in ( |2.5|) is of the form 

( 8 \ 2 

dtfsinfl 



2tt 



T / ^ 

47T JO 




+ 



2tt 



dip / d^sin^ 



i+ 



where v = const is the null congruence parallel to X + : 
(Vt») 2 = 0, (Vu, Vu) \x+= -2. 



(A.30) 



(A.31) 



The derivatives dCi /du and dC} /du of the functions appearing in (|A.21|) , (|A.22|) 
and ( |A.30[ ) are the Bondi-Sachs news functions. They determine the manifestly 
positive energy flux of the gravitational radiation. They are also the data (final data 
in the case of X + , initial data in the case of X~) for the two degrees of freedom of 
the gravitational field counted for a three-dimensional point. The three-dimensional 
points for which the degrees of freedom of the asymptotically flat gravitational field 
are counted are points of X + or X~. Finally, the last term in ( |A.30|) is the energy 
flux of the source. It will be non-vanishing at X + only if the source contains a 
massless component. 

The metric near X~ is of a form similar to ( |A.15|) but, in this case, one considers 
a converging null congruence 



(Vv) 2 = , (Vv, Vr) \ x - 



(A.32) 
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and the limit of infinite luminosity distance r is reached by going along the light 
rays towards the past. Equations similar to ( A.20|) - (|A.22 ) and ( |A.26|) — ( A.27Q define 



the Bondi mass M_(t>) and news functions 9C{" /dv, dC^ /dv at T~ . The Bondi 
mass M_(v) is the energy brought to the system by the instant v of advanced time, 
and its limit at v — > +00 is the full ADM mass: 

M_(+oo) = M . (A.33) 

The difference M — M_(v) is the energy brought to the system by the incoming 
radiation after the instant v of advanced time. Again, this is a sum of the energies 
of the incoming gravitational radiation and of an incoming flux of massless sources. 
The energy conservation law at Z~ is obtained from (|A.30| ) by changing the notation 
and the sign of the left-hand side. 

Of interest are also the limits M_(— 00) and M + (+ 00). Since we have 

dM_ 



Mn = M_ 



r°° divi 

-00) + / d^ — - , (A.34) 

J-00 dv 




M = M + (+oo) + / du — — - (A.35) 



and the integrals on the right-hand sides are, respectively, the energy brought by 
radiation and the energy carried away by radiation during the whole history, we 
conclude that M_(— 00) is the energy brought to the system by massive (time-like) 
sources coming from I~ . Similarly M + (+oo) is the remainder that goes to I + with 
time-like sources after the total emission of radiation. 

By assuming that T^ u satisfies the dominant energy condition, it has been 
proved that not only the radiation fluxes, both gravitational and matter, are 



non-negative, but also that the ADM mass and Bondi masses are always non- 
negative, which means, in particular, that the asymptotically flat system cannot 
radiate more energy than initially stored. 

The functions and A/+ appearing in eqs. ( |A.23j ) and ( |A.24j) and the func- 



tions J\f 9 and M- similarly defined at X~, are associated with radiated angular 
momentum. 

To calculate the surface integral fl2.12j ) over the boundary 
rul + : t(x) = (A.36) 
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one may first go over to the Penrose space where X~ and X + really exist as three- 
dimensional null surfaces, but one may also do this calculation directly in the phys- 
ical space-time by considering the null surfaces close to infinity and, next, going 
over to the limit of Z~ or X + respectively The contribution of X + is then of the 
form 



s 



K \I+ 



£8 16^ 



d xJ—g S 



v(x) 



+ c A»V L 



v(x) 



c) , (A.37) 



with v(x) defined in eq. ( |A.31| ), and 

= g aP V»h a p - g m V p h aP . (A.38) 

We shall carry out the calculation of ( |A.37| ) by using the Bondi-Sachs frame, 
where the metric has the form ( |A.15| ). Accordingly, the flat metric g^ v has the form 

g^dx^dx" = -du 2 - 2dudr + r 2 (d9 2 + sin 2 9d^ 2 ) , (A.39) 

and the components of h a p that do not vanish identically can be read off from 



huu 
h 

h u e 
h 

hee 

h^pi^p 

ha, n 



2GM- 



h ru 
he u 
hipu 



O 



-2Ml + O 



-2 sin 2 # A/? + O 



'2Cf 



r sin 9 



2Cf 



+ O 



h 



tp6 



1 ■ n ( ^ C2 

r sin 9 



(A.40) 
(A.41) 
(A.42) 
(A.43) 
(A.44) 

(A.45) 

(A.46) 



In the asymptotic domain near X + , eqs. ( |A.31|) are solved by 

v(x) = 2r + u, r — > 00 at fixed u. (A. 47) 

This makes it possible to carry out the integral over r in ( |A.37| ) with the aid of the 
delta-function, and the result is the limit r — > 00 of the remaining three-dimensional 
integral: 



S 



K 



1+ 



lim , 

r^oo iQ n G Jo 



2k 



dip / dflsinfl / dur 2 A»V 



r H — u 
2 



(A.48) 
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For the integrand in ( A.48|) we have 



A^Va \r + -u 



S? - K) (V" M . (A.49) 



By calculating the covariant derivatives in the metric (|A.39 ), one obtains 



-d r h, 



h, 



■In 



h 



r 6 sin V 



(A.50) 



2 2 1 

d r h uu h"uu ~t~ &r four &uh"ur ~t~ ~t~ 2 ^O^ud 



1 COS0 

H — 5- — — 77 + 

r z smu 



2-2/) u p"'u<p ; 

r z sin 



and 



Q al3 h a p — —2 ft, ur + — /lee + 9 . 2 7; ^ w 



(A.51) 



(A.52) 



It is now easy to insert the asymptotic behaviors ( A.40|) -( A.46 ) to get, for ( A.49|) : 



lim r 2 A^V u \ r + -u 

r— >oo 



: -2 GM+ - 2 - 2 (sin 6Nl) . (A.53) 

sin 6* v 7 

Upon integration over the sphere in ( |A.48|) , the contribution of J\f+,J\f_! vanishes, 
and we obtain 



5 



A" 



1+ 



1 /"27T /"7T /"OO 

/ cfy> / dflsintf / duM 

8n Jo Jo J -00 



+(u,<p,0) 



(A.54) 



A similar result, but with the opposite sign, is obtained at I . The final expression 
for Sjclg] can thus be presented in the form 



Sk = -— ( / M- 

8n \Jx+ 



M- 



duM 4 



u 



dvM_[ 



(A.55) 



where M + (u) and M_(v ) are the Bondi masses at X + and X , and the retarded and 
advanced time are normalized by the conditions ( A.18 ) and ( A.32 ). 



43 



References 

[1] D. Amati, M. Ciafaloni and G. Veneziano, Phys. Lett. 197B (1987) 81; Int. J. Mod. 
Phys. A3 (1988) 1615. 

[2] D. Amati, M. Ciafaloni and G. Veneziano, Phys. Lett. B 216 (1989) 41; 

G. Veneziano, in Superstring 89, Texas A. and M. University 1989, eds. R. Arnowitz 
et al. (World Scientific, Singapore 1990) p. 86. 

[3] G.'t Hooft, Phys. Lett. 198B (1987) 61; 

see also C. Klimcik, Phys. Lett. B208 (1988) 373. 

[4] I. Muzinich and M. Soldate, Phys. Rev. D37 (1988) 353. 

[5] E. Verlinde and H. Verlinde, Nucl. Phys. B371 (1992) 246. 

[6] A. Bellini et al., Nucl. Phys. B356 (1991) 69. 

[7] D. Amati, M. Ciafaloni and G. Veneziano, Nucl. Phys. B347 (1990) 550; Nucl. Phys. 
B403 (1993) 707. 

[8] See, e.g., R.F. Stark and T. Piran, in Proc. of the 14th Yamada Conference on 
Gravitational Collapse and Relativity, Kyoto 1986, eds. H. Sato and T. Nakamura 
(WSPC, Singapore 1986) p.249. 

[9] PD. D'Eath, Phys. Rev. D18 (1978) 990 

[10] PD. D'Eath and P.N. Payne, Phys. Rev. D46 (1992) 658, 675 and 694. 

[11] D. Amati, M. Ciafaloni and G. Veneziano, Phys. Lett. B289 (1992) 87. 

[12] L.N. Lipatov, Nucl. Phys. B365 (1991) 614. 

[13] G.W. Gibbons and S.W. Hawking, Phys. Rev. D15 (1977) 2752. 

[14] R. Penrose, in Relativity, Groups and Topology, eds. C. DeWitt and B.S. DeWitt 
(Gordon and Breach, New York 1964) 

R. M. Wald, General Relativity (The University of Chicago Press, Chicago 1984) 

[15] R. Sachs, in Relativity, Groups and Topology, eds. C. DeWitt and B.S. DeWitt (New 
York: Gordon and Breach, 1964) 

[16] E.S. Fradkin, Proc. Int. Winter School on Theoretical Physics at JINR (Dubna, 
1964); Acta Phys. Hung. XIX (1964) 175. 

[17] V. Fock, Phys. Zeit. Sowet. Un. 12 (1937) 404; 
Y. Nambu, Prog. Theor. Phys. 5 (1050) 82; 
J. Schwinger, Phys. Rev. 82 (1951) 64. 



44 



[18] R. P. Feynman, Rev. Mod. Phys. 20 (1948) 327. 

[19] G.A. Milekhin and E.S. Fradkin, Soviet Physics JETP 18 (1964) 1323. 

[20] See, e.g., B.M. Barbashov, Soviet Physics JETP 21 (1965) 402; 

V.N. Pervushin, Theor. Math. Phys. 4 (1970) 643 and references therein. 

[21] D. Kabat and M. Ortiz, Nucl. Phys B388 (1992) 570. 

[22] H. Cheng and T.T. Wu, Phys. Rev. Lett. 22 (1969) 666; 

H. Abarbanel and C. Itzykson, Phys. Rev. Lett. 23 (1969) 53. 

[23] PC. Aichelburg and R.U. Sexl, Gen. Rel. Grav. 2 (1971) 303. 

[24] See, for example: L.D. Landau and E.M. Lifshitz, The Classical Theory of Fields 
(Pergamon Press, Oxford 1975). 

[25] R. Iengo and K. Lechner, Nucl. Phys. B335 (1990) 221; 
G. Cristofano et al., Phys. Lett. B246 (1990) 45; 
M. Fabbrichesi and R. Iengo, Phys. Lett. B264 (1991) 319. 

[26] V. Ferrari, P. Pendenza and G. Veneziano, Gen. Rel. Grav. 20 (1988) 1185. 

[27] T. Dray and G. 't Hooft, Nucl. Phys. B253 (1985) 173. 

[28] C. Darwin, Proc. R. Soc. A249 (1959) 180. 

[29] N. Sanchez, Phys. Rev. D18 (1978) 1030 and 1798; 

see, also, J.A.H. Futterman, F.A. Handler and R.A. Matzner, Scattering from Black 
Holes (Cambridge University Press, Cambridge 1988). 

[30] J.B. Hartle and S.W. Hawking, Phys. Rev. D13 (1976) 2188. 

[31] E.S. Fradkin and G.A. Vilkovisky, Phys. Rev. D8 (1973) 4241. 

[32] E.S. Fradkin and G.A. Vilkovisky, 8th Intern. Conf. on General Relativity and Grav- 
itation, Abstracts of Contributed Papers (University of Waterloo, Canada 1977) pp. 
148-150. 

[33] E.S. Fradkin and G.A. Vilkovisky, Phys. Lett. B73 (1978) 209. 

[34] G.W. Gibbons, S.W. Hawking and M.J. Perry, Nucl. Phys. B138 (1978) 141. 

[35] R.M. Schoen and S.-T. Yau, Phys. Rev. Lett. 42 (1979) 547. 

[36] See R.M. Wald, ref. [14], and references therein. 



45 



Figure Captions 



Fig. 1: Tree diagrams (after cutting the scalar propagators) which are re- 
summed by our method in the case of elastic scattering. The blob represents all 
(including disconnected) trees. 

Fig. 2: Same as in fig. 1 for the case of inelastic processes. 
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